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function [cost, out 1, out 2] =cumulant_costFIRa? (par, input , pi, pll, 
p2,p22,a3,a4) ; 

% [cost, outl , out2 ] =cumulant_cost FIRa2 (par, input, pi, pll, p2 , p22, 
a3, a 4 ) ; 

% cumulant cost 

% FIR representation used 

% filter function used in both directions (non_causal) 

[np,mp]=size(par) ; 
firl=par (l:pl) ; 
. firll=par (1+pl :pl+pll) ; 
f ir2=par ( 1+pl+pll : pl+pl l+p2 ) ; 
f ir22=par (pl+pll+p2+l : mp) ; 

den=l; %FIR only 

out 1= [input ( : ; 1) -filter (f irl, den, input ( : , 2) ) -f lipud ( filter 
([0 firll ], [den], f lipud (input <:, 2))) )]; 

%/std (input (1, : ) ) ; %dlsim 

%filter 

out2= [input ( : , 2) -f ilter ( f ir2 , den, input ( : , 1) ) -f lipud ( filter 
([0 fir22 ], [den], f lipud (input < : , 1) ) ) ) ] ;%/std (input (:, 2)) ; 

%dlsim %filter 

out=[outl out2]; 

%outl=out 1/std (outl) ; % this scaling was not needed in examples ' 
in SIP98 paper 
%out2-out2/std (out2) ; 

Ld=0; % number of delays in calculating the cross-correlation 

cost3=0 ; 

cost4=0; 
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costALLl^ t) ; 
costALL2=U ; 
ol2=outl. *out2; 
cost2=mean (ol2) A 2; 

oll2=outl . *outl. *out2; 
o!22=outl . *out2 . *out2; 

if a3 ~1 

cost3=[mean(ol22) ] *2+ [mean (ol22) ]^2; 

end 

if a4 ==1 

cost4=[mean(oll2. *outl) -3*mean (outl. *2) *mean(ol2) ] *2-K . . 
(mean( (outl. ^2) ,*(out2.~2) ) -2*mean (ol2) "2- 
mean (outl . "2) *mean (out 2, A 2) J ~2+. . . 
[mean (ol22. *out2 > -3*mean (out2. A 2) *mean (ol2) ] ~2 ; 

end 

%cum4a= [cum4x (outl , outl , outl, outl ) ] A 2; 
%cum4b=[cun\4x(out2,out2,out2,out2) ] A 2; 

Cost=cost2+a3*cost3+a4*cost4; %-cum4a-curn4b; 
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